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Abstract. Around several single DAZ and DBZ white dwarfs metal-rich disks have been ob- 
served, which are mostly believed to originate from disruption of smaller rocky planetesimals. 
In some cases the material does not (only) form a dusty but gaseous disk. In the case of 
SDSSJ122859.93+104032.9 the double peaked infrared Call triplet A A 8498, 8542, 8662 A, one 
of only two emission features of the spectra, exhibits a strong red/violet asymmetry. Assuming a 
composition similar to a chondrite-like asteroid, being the most prominent type in our own solar 
system, we calculated the spectrum and vertical structure of the disk using the Tubingen non-LTE 
Accretion Disk code AcDc. Modified to simulate different non axis-symmetrical disk geometries, 
the first preliminary results are in good agreement with the observed asymmetric line profile. 
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MOTIVATION 
Metal-Rich Dust Disks 

Starting with the analysis of G 29-38 by Zuckerman and Becklin [1] in 1987, about 
twenty metal-enriched white dwarfs were discovered, which show a significant infrared 
excess in their spectra, but for which no cool companions could be found. As these 
objects have rather short sedimentation times [2] due to their high surface gravity there 
has to be an accretion source like a metal-rich dust cloud to explain their high metallicity. 
In fact, Spitzer observations confirmed the dust material [3], which is located in an 
equatorial plane [4] forming a metal-rich but hydrogen- and helium-poor [5] dust disk. 



Metal-Rich Gas Disks 

Kilic et al. [6] suggested that the hottest of those white dwarfs, with T e ff > 20 000 K, 
might sublimate their dust material and feature a gaseous disk instead. In fact, five 
white dwarfs have been found in the Sloan Digital Sky Survey (e.g. Fig. 1) by Gansicke 
et al. [7, 8, 9, priv. comm.] indicating gas disk emission features in their spectra in 
addition to the IR-spectral component of a dust disk. These gas features are mainly the 
Can A A 8498, 8542, 8662 A, although at least two of the five spectra also show a Fell 
A A 5018, 5169 A in emission. 
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FIGURE 1. Spectrum of SDSS J1228+1040 from the Sloan Digital Sky Survey. 

Origin of the Disk Material 

Theories predict that former planetary systems around the white dwarfs progenitor 
have a chance to survive the host stars' late evolutionary phases (e.g. [10]). After the 
actual formation of the white dwarf, a smaller rocky body like an asteroid or planetesimal 
might have been disturbed in its orbital motion by another surviving, larger planet and 
eventually reached the inner part of the system. Close to the central object, tidal forces 
will disrupt the asteroid providing the material to form the dust and gaseous disk. 

AIM AND METHOD 

Asymmetric Line Profiles 

Taking a closer look at the prominent Call feature of SDSS J122859. 93+104032. 9 
(hereafter SDSS J1228+1040) (Fig. 2) one recognizes the expected double-peak struc- 
ture of the lines due to the Kepler rotation of the emitting disk. But unexpectedly, the 
lines show a significant difference in the red/violet peaks line strengths. We assumed 
this asymmetry to directly originate from an actual asymmetric disk geometry. To inves- 
tigate this further we simulated a hydrodynamic disk-evolution scenario to get plausible 
geometries and later modified our spectrum- simulation code accordingly. 
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FIGURE 2. Can XX 8498, 8542, 8662 A in the spectrum of SDSS J1228+1040. Each line is splitted 
due to the Kepler rotation of the disk but shows an asymmetry in the relative strengths of the red/violet 
peaks, favoring the red wing peak in this case. 



TABLE 1. FARGO simulation parameters. 



White dwarf mass 


M W d/M 


0.77 


White dwarf radius 


# WD /km 


7700 


Disk mass 


Mlisk/g 


7 • 10 21 


Disk outer radius 


^o/%D 


136 


Simulation rim 


^sim/ 


1.5 



Hydrodynamic Simulation with FARGO 

To simulate the geometrical evolution of a gaseous disk we used the code Fast 
Advection in Rotating Gaseous Objects (FARGO) by Masset [11]. Although this hydro- 
dynamic code is suited for all kinds of sheared fluid disks, it is almost exclusively used 
in the calculation for planet formation scenarios. We calculated the temporal evolution 
of a blob of material with Gaussian density distribution (parameters given in Tab. 1). 
As a result of this, we selected two geometries, a spiral arm structure of the early simula- 
tion phase and a fully closed but out-of-balance disk from the late stage of the evolution 
(Fig. 3) to perform the modifications on the spectrum synthesis code. 



Spectral Synthesis Code AcDc 

For the calculation of the disk spectra we used the Accretion Disk code (AcDc) 
developed by Nagel et al. [12]. Starting with the assumption of a thin a-disk as described 
by Shakura and Sunyaev [13], which allows us to decouple the vertical and radial 
structure, the disk is separated into concentric annuli. For each of these rings with radius 
R the effective temperature T e ff as well as a viscosity times surface mass density wZ 
value can be derived for given white dwarf radius 7?wd, mass Mwd, accretion rate M, 





FIGURE 3. The two selected evolution phases of the FARGO simulation - color-coded is the surface 
mass density on a logarithmic scale. The left panel shows an early stage (t s { m « 40 yr) as the disks forms a 
spiral-arm like structure. The right panel shows a late (t s [ m « 125yr), but with several hundreds of orbits, 
long-lasting phase of a fully closed disk with an out-of-balance matter distribution. 
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The program then simultaneously solves the following set of equations for each of the 
rings: 

• Radiation transfer for the specific intensity / (v, jli), 

• Hydrostatic equilibrium between gravitation, gas and radiation pressure, 

• Energy conservation for viscously generated 2? mec h and radiative E m a and 

• Static rate equations ^ = for the population numbers ni of the atomic level i of 
the given model atom in the non-LTE case. 

The actual modification of the code took place in the subsequent spectral surface inte- 
gration phase of AcDc. We constructed simple geometry maps (Fig. 4) following the 
selected FARGO pictures and reset the spectral flux for those ring segments, which are 
not part of the structure to zero. AcDc then integrates the whole disk spectrum by simple 
linear interpolation between the neighboring annuli spectra. 



SYNTHETIC SPECTRA AND RESULTS 



Calculated with the set of parameters given in Tab. 2 and the model atoms (Tab. 3), the 
spectra of the two geometry maps are shown in Fig. 5. 
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FIGURE 4. Geometry maps representing the selected FARGO phases shown in Fig. 3. The two colors 
(grey shades) indicate, that the flux of the adjoint ring segment is set to zero and AcDc will compute a 
downward interpolation. For the subsequent spectrum calculations the observer is meant to look from the 
lower end of the box towards the central object. 



The first result of the analysis is, that the inner radius of the disk must be increased to at 
least 7?i>58/?wD in order to suppress broad line wings and to fit the steep line profiles as 
well as the flat continuum between the triplet components. This value is slightly larger 
than the latest results by Melis et al. [14] with R { = 40 ± 37? WD . 

The main goal of our work, matching the asymmetric line profile, also was achieved 
well for both geometries, although the relative strengths of the red/violet parts of each 
line seem to be generally estimated too low in the case of the closed disk structure. 
One has also to keep in mind, that the asymmetry strongly depends on the orientation of 
the thicker region towards the actual line of sight, so the line profile should vary with the 
rotation of the asymmetrical disk around the central object. Even though Gansicke et al. 
[9, priv. comm.] reported on a change in the line profile in the order of years for at least 
two of the known five objects, the predicted variability timescale of our model seems to 
be much shorter (order of hours) than the observations do suggest. 

TABLE 3. Model atoms and element abundances. 

Ions Hi-ii Ci-iv Oii-iv Mgi-m Sii-iv Cai-iv 

NLTE levels 11 45 32 57 40 38 

Line transitions 45 88 65 104 49 64 

Abundance (%, in mass frac.) 10" 8 4.6 65.5 13.5 15.1 1.3 



TABLE 2. AcDc integration parameters. 



Ring radii range R / 7?wd 

Disk surface mass density £/ (g/ cm 2 ) 
Disk temperature range T e ff / K 



2 - 136 
0.3 

6700 - 5600 
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FIGURE 5. Synthetic spectra for both geometries (shown in Fig. 4), in comparison with the observation. 
The line asymmetry as well as the flat continuum between the lines are fitted well for an inclination of 
i — 76°. The marked C II feature seen to the far right of the modeled spectrum is actually a yet unsplitted 
doublet line. 



On the other hand, the FARGO simulations suggest that there might also be a long-time 
variability as the disk evolves through different geometries in the time of hundreds of 
orbits. 
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